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Abstract 

Dynamics  of  ideal  fluid  with  free  surface  can  be  effectively  solved  by 
perturbing  the  Hamiltonian  in  weak  nonlinearity  limit.  However  it  is 
shown  that  perturbation  theory,  which  includes  third  and  fourth  order 
terms  in  the  Hamiltonian,  results  in  the  ill-posed  equations  because  of 
short  wave  instability.  To  fix  that  problem  we  introduce  the  canonical 
Hamiltonian  transform  from  original  physical  variables  to  new  variables 
for  which  instability  is  absent.  We  found  the  choice  of  such  transform  is 
unique. 

1  Introduction 

The  Euler  equations  describing  dynamics  of  ideal  fluid  with  free  surface  is  a 
Hamiltonian  system,  which  is  especially  simple  if  the  fluid  motion  is  potential, 
u  =  V$,  where  v  is  the  fluid’s  velocity  and  $  is  the  velocity  potential.  In  this 
case  [1,  2]  the  Euler  equations  can  be  presented  in  the  form: 

dv  =  SH  _  _6H 

dt  <5$’  dt  ~  Stj  '  ^ 

Here  z  =  J?(r)  is  the  shape  of  surface,  z  is  vertical  coordinate  and  r  =  (x.  y)  are 
horizontal  coordinates,  $  =  is  the  velocity  potential  on  the  surface.  The 
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Hamiltonian  H  coincides  with  the  total  (potential  and  kinetic)  energy  of  fluid. 
The  Hamiltonian  cannot  be  expressed  in  a  closed  form  as  a  function  of  surface 
variables  t},  $,  but  it  can  be  presented  by  the  infinite  series  in  powers  of  surface 
steepness  |  Vr/|: 

H  =  Ho  +  Hi  +  H2  +  —  (2) 

Here  H0,  Hu  H%  are  quadratic,  cubic  and  quartic  terms,  respectively.  Equa¬ 
tions  (1),  (2)  are  widely  used  used  now  for  numerical  simulation  of  the  fluid 
dynamics  [7,  8,  9,  10,  12,  13,  14,  15],  These  simulations  are  performing  by  the 
use  of  the  spectral  code,  at  the  moment  a  typical  grid  is  512  x  512  harmonics. 

Canonical  variables  are  used  also  for  analytical  study  of  the  surface  dynamics 
in  the  limit  of  small  steepness.  It  was  show  [4,  5,  6]  that  the  simplest  truncation 
of  the  series  (2),  namely 

H  =  H0  +  H1)  (3) 

leads  to  completely  integrable  model  -  complex  Hopf  equation.  In  framework  of 
this  approach  one  can  develop  the  self-consistent  theory  of  singularity  formation 
in  absence  of  gravity  and  capillarity  for  two  dimensions  (one  vertical  coordinate 
z  and  one  horizontal  coordinate  x). 

However,  canonical  variables  rj,  $  has  a  weak  point,  which  becomes  clear,  if 
we  concentrate  our  attention  on  the  complex  Hopf  equation, 

dQ  _  1  /c?@\2 

~dt  ~  ~2\dx)  ’  (4) 

which  comes  from  Eqs.  (1),(3).  Here 

*  =  &(©)  (5) 

and  ©  is  analytic  function  of  complex  variable  £  in  a  strip  ~h  <  Im(x)  <  0. 
The  weak  point  is  that  Eq.  (4)  is  ill-posed.  A  general  complex  solution  of  this 
equation  is  unstable  with  respect  to  grow  of  small  short-wave  perturbations.  The 
same  statement  is  correct  with  respect  to  more  exact  fourth  order  Hamiltonian 

H  =  Ho  +  Hi  4-  (6) 

which  is  used  in  most  numerical  experiments.  These  experiments  are  easy  be¬ 
comes  unstable:  to  arrest  instability  one  should  include  into  equations  strong 
artificial  damping  at  high  wave  numbers.  Even  in  presence  of  such  damping 
one  can  simulate  only  waves  of  a  relatively  small  steepness  (not  more  that 
0.15)  [reference??????????] . 

In  this  Article  we  show  that  these  difficulties  can  be  fixed  by  a  proper  canon¬ 
ical  transformation  to  another  canonical  variables.  It  is  remarkable,  but  the 
property  of  nonlinear  wave  equation  to  be  well-  or  ill-posed  is  not  invariant  with 
respect  to  choice  of  the  variables . 

In  the  present  Article  we  demonstrate  that  there  are  uniQue  canonical  vari¬ 
ables  such  that  the  Eqs.  (1),  (6)  are  well-posed.  We  call  these  variables  “optimal 
canonical  variables”.  We  can  formulate  a  conjecture  that  the  optimal  canonical 
variables  exist  and  are  unique  in  all  orders  of  nonlinearity. 
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2  Basic  equations  and  Hamiltonian  formalizm 

Consider  the  dynamics  incompressible  ideal  fluid  with  free  surface  and  constant 
depth.  Fluid  occupies  the  region 

- h<z<p(r ),  r  =  (a;,y),  (7) 

where  (a;,  y)  are  the  horizontal  coordinates  and  z  is  vertical  coordinate. 

Viscosity  is  assumed  to  be  absent  and  the  fluid’s  velocity  v  is  potential  one: 

v  =  V$,  (8) 

where  $  is  the  velocity  potential.  Incompressibility  condition, 


V  *  v  =  0,  (9) 

results  in  the  Laplace  Eq. 

=  0.  (10) 

The  potential  satisfies  also  the  Bernoulli  equation: 

$t  +  ^(V$)2+P  +  02  =  O,  (11) 

where  p  is  the  pressure,  g  is  the  acceleration  of  gravity,  and  we  set  density  of 
fluid  to  unity. 

There  are  two  types  of  boundary  conditions  at  free  surface  for  Eqs.  (10),  (11). 
First  is  the  kinematic  boundary  condition 

J  =  (*,  -  Vr,  •  V$)  |  =  *Wl  +  (V*?)2,  (12) 

where  vn  =  n  •  V#  is  the  normal  component  of  fluid’s  velocity  at  free  surface, 
and  n  =  (-V77, 1)  [l  +  (Vf?)2]_1//2  is  the  interface  normal  vector. 

Second  is  the  dynamic  boundary  condition  at  free  surface 


■ 


Vt? 


\/i+W’ 


(13) 


where  a  is  the  surface  tension  coefficient  which  determines  the  jump  of  pressure 
at  free  surface  from  zero  value  out  of  the  fluid  to  p\z_  value  according  to  Eq. 
(13). 

Boundary  condition  at  the  bottom  is 


ft* 


(14) 


Eqs.  (10)  -  (14)  form  a  closed  set  of  equations  to  determine  the  dynamics 
of  free  surface. 
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The  total  energy,  H,  of  the  fluid  consists  of  kinetic  energy,  T,  and  potential 
energy,  U : 

H  =  T  +  U,  (15) 

T=lJdrJ\v$)2dz,  (16) 

u=\9 1 1,2  dr+ °  /  bfcw?  -  i]dr'  w 

It  is  a  convenient  to  introduce  the  value  of  velocity  potential  at  interface 
with  the  boundary  conditions 


$|^  =  *(r,f).  (18) 

It  was  shown  in  Ref.  [2]  that  free  surface  problem  (10)  —  (14)  can  be  written  in 
the  Hamiltonian  form  (1),  with  the  Hamiltonian  H  defined  in  (15). 

Fourier  transform, 

$k  =  -^  J  exp(-ik  •  r)$(r)dr,  (19) 

is  the  canonical  transform  which  conserves  the  Hamiltonian  structure  and  Eqs. 
(1)  take  the  following  form: 


drjk  _  SJI  89 k 
dt  69- k  ’  dt 


5H 

*?-k’ 


Vk 


j?_k,  K  =  *-k. 


(20) 


3  Weak  nonlinearity 

If  typical  slope  of  free  surface  is  small,  |  Vt?|  <  1,  the  Hamiltonian  H  is  can  be 
series  expanded  (see  Eq.  (2))  in  powers  of  stepness  |Vt/|  which  gives  [2,  3]: 

H0  =  -  /  {^|^k|2  +  Bfc|r/k|2  jdk, 

Ak  =  k  tanh(fch),  Bk  =  9  4*  crfc2,  k  =  |k|,  (21) 

Hl  =  ~bJ  W(ki  +k2+  k3)dkidk2dk3  (22) 

~  2(2ir)2  f  [£'ki),ka1k3,k4$ki^,k2  -  ^(kj  •  k2)(k3  •  k4)r/kl7?k2] 

X-VksVkAki  +  k2  +  k3  +  k4)dkidk2dk3dk4,  (23) 
where  matrix  elements  are  given  by 

LW t*2  =  ~kl  ‘ ka  “  AiAq’ 

•C'k1),k2,k3,k4  =  \AA2  (Ai+3  +  A2+3  +  A1+4  +  A2+4) 

-\{k\A2  +  %AA  Aj=Akj.  (24) 
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The  corresponding  dynamical  equations  follow  from  (1),  (6),  (21),  (22),  (23) : 


™  =  -m  +  <rAV  +  \ [(A*)2  -  (v$)2]  -  (A*)A[i,(A#)] 

-(A*) (A*)i»  -  |  V  •  [Vtj(Vtj  •  V*j)] ,  (25) 

^  =  i#  -  V  •  [(V*)»j]  -  .+  A{r)A[riAv] } 

+iA[7?2i$]  +  ii[»?2A$]1  (26) 


where  A  is  the  linear  integral  operator  which  corresponds  to  multiplication  on 
k  tanh(fch)  in  Fourier  space.  For  two  dimensional  flow,  $  (a,  y)  =  $(ar),  t)(x,y)  = 
Tj(x),  this  operator  is  given  by 


*/(*)  =  ThP-v. 


/+oo 

; 

-oo  1 


/(*') 


sinh  [(s'  -  x)Tr/(2h)] 


dx* 


(27) 

(28) 


where  P.V.  means  Cauchy  principal  value  of  integral.  In  the^  limiting  case  of 
infinitely  deep  water,  h  ->  oo,  operator  A  turns  into  operator  k 

lim  A  —  k  (29) 

h-+oo 

which  corresponds  to  multiplication  on  |k|  in  Fourier  space  while  operator 
R  for  two-dimentsional  flow  turns  into  the  Hilbert  transform: 

lim  R  =  H,  Hf(x)  =  -P.V.  f+°°  4^-dx'.  (30) 

h — hoo  It  y_00  T  -  X 


H  can  be  also  interpreted  as  a  Fourier  transform  of  -isign(fc). 

If  one  neglects  gravity  and  surface  tension,  g  =  0,  a  =  0,  than  Eqs.  (1),  (2) 
at  leading  order  over  small  parameter  \Vr)\  result  in[5,  4,  6] 

g  =  A*,  (31) 

^=l[(i$)2_(V$)2].  (32) 

Remarkable  feature  of  Eqs.  (31),  (32)  is  that  the  second  Eq.  (32)  does 
not  depend  on  rj  thus  one  can  first  solve  (32)  and  then  find  rj  from  Eq.  (31). 
Substitution  0  =  $  +  iRV  into  Eq.  (32)  results  in  the  complex  Hopf  Eq.  (4) 
for  two-dimensional  flow  [6]  which  is  completely  integrable. 

Both  Eqs.  (32)  and  (4)  are  ill-posed  because  they  have  short  wavelength 
instability  which  determines  as  follows.  E.g.  we  can  analyze  Eq.  (32).  Take  4* 
in  the  form 

*  =  $0  +  (*ieik°-r+,/f  +  c.c.) ,  (33) 
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where  ^o(r,  t)  is  a  solution  of  Eq.  (32),  is  the  amplitude  of  small  pertur¬ 
bation,  and  c.c.  means  complex  conjugation.  Then  in  the  limit  |ko|  ->  oo 
evolve  very  slow  in  space  compare  to  eik°'r+l/t  and  we  get  the  dispersion  relation 
for  small  perturbations: 


v  =  AkoA^0  -  -  V^0  (34) 

which  describes  instability  for  Re(v)  =  AkoA^0  >  0.  For  general  initial  condi¬ 
tion  such  instability  region  always  exists.  Instability  growth  rate,  Re{v)  grows 
as  |k0|  increases. 

4  IU-posedness  of  the  fourth-order  Hamiltonian 

Consider  now  a  more  general  case  of  nonzero  g  and  a  and  take  into  account  all 
terms  in  the  Hamiltonian  up  to  forth  order,  i.e.  consider  full  Eqs.  (25).  Similar 
to  previous  section  we  linearize  Eqs.  (25)  using  ansatz 

r)  =  m  +  (rh^°-I+vt+c.c), 

$  =  $0  +  ($1e<kor+^  +  C.C.),  (35) 

were  r?o(r,t),  ^o(r?^)  are  solutions  of  (25),  and  get  for  |k0|  -4  oo  : 

V  ~  2  +  2  +  +  “  k2)(-l  +  Akor}o)r}o] 

x[-2«/+{3(V%)2-2}fcV])  , 

M  =  -iko  •  Vtf0  +  [Ako  +  (*S  -  4„)»?o]i«o.  (36) 

Instability  growth  rate  take  more  compact  form  for  infinite  depth  fluid  as 
limA-K»  Ak  =  k: 

u=  t  +  +  2ko[-2g+ {Z(Vt]o)2 -2}k2a]j  '  , 

fi  =  — iko  *  V^o  +  kok^o-  (37) 

For  general  initial  condition  with  arbitrary  depth  h  the  instability  region 
Re{v)  =  AkoA^o  >  0  always  exists.  We  conclude  that  full  fourth  order  Hamil¬ 
tonian  does  not  prevent  short  wavelength  instability  and  Eqs.  (25)  are  ill-posed. 

Hl-posedness  of  Eqs.  (20),  (21)  -  (23)  (or,  equivalently,  Eqs.  (25))  makes 
them  difficult  for  numerical  analysis.  There  a  few  ways  to  cope  with  that  prob¬ 
lem.  One  way  is  to  introduce  artificial  damping  for  short  wavelengths,  i.e.  to 
replace  Eqs.  (20)  by 

m  =  +  7l(fc)r?k’  ~df  =  + 72(*)$k’  (38) 
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where  functions  7 i(fc),  72  (k)  are  zero  for  small  and  intermediate  values  of  k  but 
they  tend  to  -00  for  k  00.  Second  way  is  to  introduce  finite  viscosity  of 
the  fluid.  However  in  that  case  the  Hamiltonian  does  not  conserve  and  we  can 
not  use  the  Hamiltonian  formalism.  In  this  Paper  we  use  third  way  which  is  to 
introduce  new  canonical  transform  of  variables  77,  to  remove  short  wavelength 
instability.  Advantage  of  this  methods  is  that,  in  contrast  to  first  way,  we  do  not 
introduce  any  artificial  damping,  and,  in  contrast  to  second  way,  we  preserve 
the  Hamiltonian  formulation  of  free  surface  dynamics  problem. 


5  Canonical  transform 

Canonical  transform  from  variables  77  to  new  variables  R,  £  is  determined  by 
the  generating  function  5  : 

S  =  J  +  —  J  +^2  +k3)dk1dk2rfk3 

"^4(27r)2  f  ^i»ka»^ik4^i7?ka*?^^k4^(ki  +  k2 +ka +k4)dkidk2dk3dk4,  (39) 


SS  1  f 

*k  =  Srj^  ~Rk  +  4'i r  J  AlRk^S^  +k2  -  k)dki<flc2 

+  4(2  w)2  /  ^ki,k2,k3,-k-Rki»?k2%s<5(ki  +  k2  +  k3  -  k)dkidk2dk3,  (40) 
.  SS  1  f 

&  —  §R_k  ~Ttk  +  sit  J  AkTfinVkiSiki  +  k2  -  kjdkjtfl^ 

+4(^)2  /  v'-k,k2,ks,k4»7k2»?k3»?k4^(k2  +  k3  +  k4  -  k)dk2dk3dk4,  (41) 

where  is  41x6  symmetric  function  of  k2)k3,  k4.  This  is  the  most  gen¬ 

eral  form  of  canonical  transform  up  to  terms  of  fourth  order.  The  only  condition 
which  we  use  here  is  that  S  is  chosen  to  be  linear  functional  R  to  preserve  the 
quadratic  dependence  of  the  Hamiltonian  on  canonical  momentum  R. 

i)  can  be  found  from  Eq.  (41)  as  the  functional  of  £  by  iterations  (here 
and  below  we  take  into  account  only  corrections  up  to  the  fourth  order  in  the 
Hamiltonian): 


%  -  &  “  ^  J  MuMh  +  k2  -  k)dkidk2  +  — ^ 

X  J  [■^*‘^1+2  —  2^r-k,kl,k2,ks]£ki£k2£k3<S(ki  +  k2  +  k3  —  k)dkidk2dk3,  (42) 

Eqs.  (40),  (42)  give: 


*k  =  Rk  + 


4^  f  AiR*itk,S(k! 


+  k2  —  k)dkidk2  + 
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J  [~  ^1^2+3  +6Fk1>k2,k3,-k]^ki^k2^k3 

x<5(ki  ■+•  k2  +  k3  —  k)dkidk2dk3,  (43) 

Using  Eqs.  (21),  (22),  (23),  (40),  (42)  we  get: 

Ho  =  \J  {A*|i?k|2  +  Bfc|&|2}dk,  (44) 

Hi  =  /[- (ki  •  k2)-Rk,  l?k2  —  —  (-4i  B\  +  A2B2  +  ^3  Hs)^ki  <^k2  j 

x£ks^(ki  4-  k2  +  k3)c/kieflc2c?k3,  (45) 

H2  =  J  {(^1  '  ^2)(-4l+2  -M  —  M)  -  fcl^2  -  k%Ai 

+  ^.41A2[2l)+3  +  j42+3  +  Ai+4  +  A2+4]  +  3[j4lVk2tk3,ki,ki 

+J42Vrki,ks,k4,k3]  ]-Rki#k3fks£k4<J(ki  +  k2  +  k3  +  k4)dkidk2dk3dk4 

+  g^2  J  {- ff(ki  -k2)(k3  -k4)  +  ^?+2Hi+2  +  A3B3Ai+2 

-2Bi  Vkj  tk2  ,k3,k4  )  £ki  £k3£k$£k4 
x  S  (ki  +  k2  +  k3  +  k4  )dki  dk2  dk3dk4 ,  (46) 

Canonical  transform  conserves  the  Hamiltonian  structure  so  the  dynamical 
equations  in  new  variables  R,  £  are  given  by: 

d$  _  5H  8R  =  _SH  ,47 s 

8t~  SR'  dt  SC 

6  From  complex  to  real  Hopf  equation 

We  chose  the  cubic  term  of  the  generation  function  S  in  such  a  way  to  remove 
linear  instability  at  leading  order.  Similar  to  Eqs.  (31),  (32),  we  get  from  Eqs. 
(44),  (45),  (47)  at  leading  order  of  small  parameter  |V£| : 

%  =  AR,  (48) 

ot 

™  =  -i(VH)2.  (49) 

at  2'  ’ 

Thus  instead  of  the  complex  Hopf  Eq.  (4)  (or  Eq.  (32))  we  got  real  Burgers  Eq. 
(49)  for  new  canonical  variable  R.  It  is  important  that  the  real  Burgers  Eq.  is 
well-posed. 

Additional  advantage  of  (49)  is  that  it  can  be  integrated  by  the  the  method 
of  characteristic  not  only  in  two  dimensions  as  Eq.  (4)  but  for  three  dimensional 
flow  also. 
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7  Removal  of  instability  from  fourth  order  term 

Next  step  is  to  remove  instability  from  the  fourth  order  terms  in  the  Hamiltonian 
(46)  by  proper  choice  of  matrix  element  V.  We  can  take  Vk1,k2,k3,k4  in  fche 
following  form: 

Vki,k2,k3,k4  =  aikl  +  a2J4i(A2+3  +  ^2+4  +  ^3+4),  (50) 

where  ot 2  axe  the  real  constants.  . 

The  dynamical  equations,  as  follows  from  (44),  (45),  (46),  (47),  (50)  are 

™  =  -m-  '-(vRf  +  + \m  +  iei(v.R)2 

•  VAR  -  ^(1  -  3«i)£ {AR)AR  -  ( J  +  3<*2)  {AR)A(ZAR) 
_^i[(ii?)2]  -  |v-  [V£(V£-V|)]  -  \tA2(Be) 

+^a2(AB()Ae  -  + ^AB(ae),  m 

g  =  AR  -  V  •  [(VE)S]  +  \v  •  [(VH)ie2]  -  Jv  •  A{eVR) 

*  -iv  .  (e2AVH)  + 

+  (\  +3a2)A[a(0*)]  +la2A[(Ae)(AR)l  (52) 

where  B  =  g  —  <jA. 

Using  ansatz 

I  =  lo  +  (4ieiko-r+,'<  +  C.C.), 

R  =  Rc  +  (RieiUT+’'t  +  c.c.) ,  (53) 

one  can  linearize  Eqs.  (51),  (51)  on  a  background  of  solution  £0(r,  t),  Bo(r,t)  a, 
and  get  for  |ko|  -4  oo  : 

„=a+^+/13)1/2,  (54) 

fix  =  — iko  •  ViZo  +  2^fco^°^°  * 

-i  [2(3ai  -  l)kl  +  (1  +  12a2)A£0]£0ARo, 

/i2  =  »ko  •  VRo  +  ^-Afco&ko  •  VJlo 
+i  [2(3ai  -  1  )kl  +  (1  +  36a2)40]Mi?o, 
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^3  g  |2(B*;0  [4Afc0$o  —  4  4-  3A*0(4a2  —  l)£o  4-  Oaifco^o] 

+6ftgcr(V€ o)2)  +^0(6a2  -  1)5*0^} {4tgf0 

+(1  +  12«2)^0^o  +  (4  -f  6aifco£oM*o  +  (6A*0a2  —  &o)A£q  j.  (55) 

To  avoid  instability  it  is  necessary  to  have  purely  imaginary  v.  Necessary  con¬ 
dition  for  that  is  that  the  expression  under  square  root,  /x2  +  /x3,  in  Eq.  (54) 
should  have  zero  imaginary  part.  It  follows  form  Eqs.  (54)  that  /x3  is  always 
real,  and  /x2  is  real  provided 

2(3cki  —  1)&q  4-  (1  4-  36a2)^.fe0  =  0.  (56) 

Second  condition  to  have  purely  imaginary  v  is  Re{\x  1)  =  0,  which  gives  the 
second  condition 


2(3aq  —  1  )&o  +  (1  +  12a2)A  \0  =  0.  (57) 

It  follows  from  Eqs.  (56),  (57)  that  in  the  limit  fco  -4  00  (remember  that  Ak0  — > 
ko  in  that  limit)  the  system  (44),  (45),  (46),  (47),  (50)  is  well-posed  provided 

<*1  =  1/6,  a2  =  0,  (58) 

which  gives  from  (46),  (50)  the  well-posed  fourth-order  Hamiltonian: 

=  8(27r)2  J  ■®C2)(^1+2  —  A2)  —  2 (&i  A2  +  h\Ax) 

[^+3  +  ^2+3  +  ^4i+4  4-  ^2+4]  | 

x*ki*ka&.6c4*(ki  +  k2  +k3  +k4)dkidk2dk3dk4 

+  8(27r)2  f  {  <T^S‘1  *  ^(ka  ‘  +  4^1+2^142  4-  fc3H3J4i_|_2 

"3Sl^}&l^^,(kl  4“  k2  +  k3  +  k4)dk1dk2dk3dk4,  (59) 

This  Hamiltonian  is  well-posed  for  any  g ,  a  (including  case  g  =  a  =  0). 

To  find  dynamics  of  free  surface,  one  can  solve  Eqs.  for  R ,  £  using  Eqs. 
(47),  (44),  (45),  (59).  This  is  the  main  result  of  this  Article.  To  recover  physical 
variables  tf,  1 ?  from  given  R,  £  one  can  use  Eqs.  (42),  (43),  (50),  (58). 

As  follows  from  Eq.  (58),  new  canonical  variables  £,  R  are  uniquely  deter¬ 
mined  from  the  condition  of  well-posedness  of  dynamical  Eqs.  in  new  variables 
up  to  the  fourth  order  in  the  Hamiltonian.  We  refer  to  these  variables  as  optimal 
canonical  variables.  For  some  extent  similar  results  were  obtained  by  Dyachenko 
and  Shamin  [18]  for  particular  case  of  two-dimensional  flow.  We  conjecture  that 
the  optimal  canonical  variables  exist  and  are  unique  in  all  orders  of  nonlinearity. 


8  Special  cases 

There  are  a  number  of  important  special  cases  of  optimal  canonical  variables. 
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8-1  Deep  water  limit 

For  h  -¥  oo  Eqs.  (44),  (45),  (59)  take  the  form 

Ho  =  \j  {fc|f*k|2  +  £fc|&|2}dk,  (60) 

Hi  =  JJj:  J  [  —  (ki  •  k2)iZki-Rk3  —  g(fclBi  +  +  fc3-B3)fki£k2] 

x^k,J(ki  +  k2  +  k3)dkidk2dk3,  (61) 
H a  =  /  {(ki  •  k2)(|k1  +  k2|  -  fa  -  fc)  -  i(fc2fc2  +  k2k x) 

+  -kifc2[|ki  +k3|  4-  |k2  +k3|  +  |ki  +  +  |k2  +k4|]| 

X-R^-Rk^kslk,  J(ki  +  k2  +  k3  +  k4)dk1<flc2dk3dk4 
+3(2^2 /{- ^(ki  •  k2)(k3  •  k4)  +  j|ki  +  k2  |2J3i+2  +  *3.83^1  4-  k2| 

-  ^  Bi  fc?  }  £kl  fk2  £k3  Ik,  5(ki  +  k2  +  k3  +  k4)dkidk2dk3dk4.  (62) 

8.1.1  Zero  gravity  and  capillarity  g  =  a  =  0 

Ho  =  \J  fc|i?k|2dk,  (63) 

Hi  =  J (kj  •  k2)fiki f?k2^k3 ^ (kj  +  k2  +  k3)«fkidk2dk3,  (64) 

H 2  =  /  {(k i '  k2)(lki  +  k2|  -  *i  -  fa)  -  §(*fe  +  fc22*i) 

+-fcik2[|ki  4-  k3|  4-  |k2  +  k3|  4-  |ki  4-  k4|  4-  |k2  4- 1^1]  | 

x  Rk ,  Rki  ^k3  Ik,  <5(ki  +  k2  +  k3  +  k4)dk1dk2dk3dk4.  (65) 

8.1.2  Zero  gravity,  g  =  0,  and  nonzero  capillarity  a  jtO 

Ho  =  \j  {fc|f?k|2  +  <rk2|&|2}dk,  (66) 

Hl=hf  [_(kl‘  k2)i?k‘fi^  "  ^  +  *f  )&^k2] 

x^k3^(ki  4-  k2  4-  k3)dk1dk2(ik3,  (67) 

Hi  =  f  {(kx  •  ka)(|k1  +  k2|  -ki-h)-  \{k\h  +  k\ki) 

+-fcifc2[|ki  +k3|  +  |k2  +k3|  4-  |k4  +k4|  4-  |k2  +  k4|]| 
x  jRk,  iik2  ?k3  <5(ki  +  k2  +  k3  4-  k4)dkidk2dk3dk4 
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+  8(ir)2  /  {  —  (ki  •  k2>(k3  •  k4)  +  ^Iki  +  lc2|4  +  fc||ki  +  k2| 

~3^i}&iCk3fk3Ck4<J(ki  +  k2  +  k3  +  k4  )dki  dk2  dk3  dk4 .  (68) 

8.1.3  Nonzero  gravity,  g  /  0,  and  zero  capillarity  a  =  0 

Ho  =  \f  {fc|-Rk|2  +  </|&|2}dk,  Bk=g  +  ak2,  (69) 

^1  =  4irf\~^Cl  k2)#ki#ka  ~  +  fag  +  fc3ff)£ki£k2] 

x^ka^kx  +  k2  +  k3)dkidk2dk3,  (70) 

^2  =  8(2 xr)2  J  {^kl  ~  fa  ~  fa)  ~  ^{klfa  +  k2k i) 

+\fafa  [|kx  +  k3|  +  |k2  +k3|  +  |kx  +  k4|  +  |k2  +  k4|] } 
xi?kli2k3Ckslk4<i(ki  +  k2  +  k3  +  k4)dkxdk2dk3dk4 

~*~8(27r)2  y  {  —  ’  ^2)(^3  ‘^£4)  +  +k2|2p  +  fc3(j|ki  +k2| 

-^fcij&j&a&s&^ki  +  k2  +  k3  +  k4)dkidk2dk3dk4.  (71) 

8.2  Shallow  water  limit 

Shallow  water  limit  corresponds  to  kh  ->•  0.  Li  that  limit  Ak  -x  k2h.  Eqs. 
(44),  (45),  (59)  take  the  form 

Ho  =  lf  {*2h|J?k|2  +  S*|&|2}dk,  (72) 

^1  =  4nf\.~^1  k^-Rki-Rka  -  g  (fcl-Bi  +  fcfi?2  +  fc3B3)|kifka] 

x?k3<i(ki  +  k2  +  k3)dkxdk2dk3,  (73) 

H 2  =  s(Ly  /  (2(kl ' ka)2  -  \^A2  + 

+— fcifcl  [|ki  +  k3|2  +  |k2  +  k3|2  +  |ki  +  k4|2  +  |k2  +  k4|2]| 

x  Rki  -Rka  £k3  ^k4  <£(kl  +  k2  +  1^3  +  k4)dJcidlc2<ik3dk4 

+  8(2xr)2  J  {  _  <T^kl  ’  ^2X^3  ‘  ^4)  +  ‘4_l^i  +  fa\4Bi+2 

+k3B3\k1+k2\2-^B1k2} 

x^k,^k2^k3^k4^(ki  +  k2  +  k3  +  kj  )dki  dk2dk3  dk4 .  (74) 

However  conditions  (56),  (57)  can  not  be  simultaneously  satisfied  for  shallow 
water.  Thus  shallow  water  problem  remains  ill-posed  even  for  optimal  variables 
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because  of  short  wavelength  instability  coming  from  the  fourth  order  term  in  the 
Hamiltonian  although  instability  from  third-order  term  is  removed  by  canonical 
transform  (40),  (41).  Actually  there  is  no  big  surprise  in  that  because  we  need 
to  choose  carefully  the  order  of  taking  the  limits  kh  -¥  0  and  k0h  -►  oo,  where 
k  is  the  typical  wavevector  of  surface  motion  and  ko  is  a  wavevector  of  short 
scale  perturbation.  It  means  that  we  need  first  to  solve  dynamics  of  water  with 
finite  depth  h  which  is  well-posed  problem  in  optimal  canonical  variables,  and 
only  after  that  we  should  take  limit  kh  -¥  0. 

9  Conclusion 

In  conclusion,  we  found  optimal  canonical  variables  for  which  the  water  wave 
problem  is  well-posed  in  the  approximation  which  keeps  terms  up  to  fourth 
order  in  the  Hamiltonian.  The  important  question  remain  open  if  it  is  possible 
to  make  water  wave  equations  well-posed  by  proper  choice  of  canonical  transform 
for  higher-order  corrections  (fifth  and  higher  order).  We  conjecture  that  such 
optimal  canonical  variables  exist  and  are  unique  in  all  orders  of  nonlinearity. 
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